Two-temperature coronal flow above a thin disk 
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ABSTRACT 

We extended the disk corona model (Meyer & Meyer-Hofmeister 1994; Meyer, Liu, & Meyer- 
Hofmeister 2000a) to the inner region of galactic nuclei by including different temperatures in ions 
and electrons as well as Compton cooling. We found that the mass evaporation rate and hence 
the fraction of accretion energy released in the corona depend strongly on the rate of incoming 
mass flow from outer edge of the disk, a larger rate leading to more Compton cooling, less efficient 
evaporation and a weaker corona. We also found a strong dependence on the viscosity, higher 
viscosity leading to an enhanced mass flow in the corona and therefore more evaporation of gas 
from the disk below. If we take accretion rates in units of the Eddington rate our results become 
independent on the mass of the central black hole. The model predicts weaker contributions to the 
hard X-rays for objects with higher accretion rate like narrow-line Seyfert 1 galaxies (NLSls), in 
agreement with observations. For luminous active galactic nuclei (AGN) strong Compton cooling 
in the innermost corona is so efficient that a large amount of additional heating is required to 
maintain the corona above the thin disk. 

Subject headings: accretion, accretion disks-galaxies: nuclei-X-rays:gaIaxies-X-rays:binaries 



1. Introduction 

It is commonly thought that AGN are com- 
posed of a central supermassive black hole and a 
thin equatorial accretion disk. Spectral features 
of AGN, the big blue bump, the soft X-ray excess 
and hard X-ray tails, and the 6.4 keV fluorescent 
iron lines show strong observational evidence for 
hot gas in the neighborhood of the inner accre- 



tion disk (e.g., Mushotzky, Done, & Pounds 1993; 
Zdziarski 1999). Soft disk photons are Compton 
up-scattered by the optically thin hot plasma and 
cause the power-law hard X-ray component in the 
spectrum. The relative strength of the hard X- 
ray component with respect to the UV-soft X-ray 
component, leading to the classification as soft or 
hard spectrum, is ascribed to the dominance of ei- 
ther the thermal emission from the accretion disk 
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or of the thermal Comptonization in a hot plasma. 
The iron lines are explained by reflection of the 
hard X-rays from the underlying disk (e.g., Fabian 
et al. 2000; Reeves et al. 2001; Gondoin et al. 
2001). 

Theoretical simulations successfully reproduce 
the observed spectra by assuming a certain frac- 
tion of gravitational energy released in the hot 
gas (e.g. Haardt & Maraschi 1991; Kawaguchi, 
Shimura, & Mineshige 2001) and/or a certain spa- 
tial distribution of the hot gas (e.g. Merloni & 
Fabian 2001a). However, without understanding 
the coexistence of disk and corona, the previous 
models require additional free parameters besides 
the mass of the black hole and the accretion rate 
for modeling the observed spectra. 

On the other hand, mass evaporation from the 
underlying cool disk to the hot corona provides 
a mechanism for the change of the accretion disk 
in the outer region to the coronal flow/ AD AF in 
the inner region. Such a model (Meyer & Meyer- 
Hofmeister 1994, Meyer et al. 2000a) describes 
the physics of the process how mass and energy 
exchange and how the hot corona exists above 
the cool disk. Further the model allows to ex- 
plain self-consistently whether and where the ac- 
cretion changes from the thin disk to the hot coro- 
nal flow/ AD AF in dependence on the given ac- 
cretion rate. The model explains the soft/hard 
spectral transitions in black hole X-ray binaries 
(Meyer, Liu, & Meyer-Hofmeister 2000b) and the 
truncation of the thin disk in soft X-ray transients 
and some galactic nuclei ( Liu et al. 1999; Liu 
& Meyer-Hofmeister 2001). It provides a scenario 
for disk evolution and outburst cycles of soft X- 
ray transients (Meyer-Hofmeister & Meyer 1999). 
In these applications of the evaporation model the 
inner disk is truncated due to a relatively low ac- 
cretion rate. 

If we apply the disk corona model to AGN, the 

inferred high mass accretion rate in these luminous 
systems makes it possible for the disk to extend 
inward to the last stable orbit since the high ac- 
cretion rate tends to prohibit the disk from being 
depleted. Then the hard X-ray radiation origi- 
nates from re-processing in the hot corona. Such 
an inner corona has very high temperature and 
low density, and hence electrons and ions are no 
longer well coupled if all the heating goes to ions 
and only Coulomb collisions transfer energy from 



ions to electrons. Photons from the thin disk are 
Compton up-scattered in the coronal region close 
to the central black hole. This process influences 
the coronal structure. Our former applications 
of the evaporation model adequate for equal ion 
and electron temperature (one-temperature coro- 
nae) would not be a good approximation anymore 
for the innermost regions of the black hole accre- 
tion system. A two-temperature treatment has to 
be developed. 

Recently a semi-analytical approach based 
on the same physics is carried out in detail by 
Rozanska & Czerny (2000). In their investigation 
the corona is regarded as a layer with constant 
pressure, which means that the frictional heating 
is higher than that in the real corona where the 
pressure scale height is much smaller as shown by 
our numerical computations (Meyer et al. 2000a). 
This approximation and further assumptions in 
their work lead to quite a difference in the de- 
tailed results. A recent investigation by Spruit 
& Deufel (2001) concerns the heating of the in- 
nermost disk region by hot ions and the resultant 
evaporation of matter from the disk. 

In the present study we took into account the 
decoupling of electrons and ions and the Compton 
cooling, extending the previous disk evaporation 
model (Meyer et al. 2000a) to small distances 
from the black hole. We calculated the vertical 
structure of the corona numerically, attempting to 
describe the hot corona self-consistently in equi- 
librium with the cool disk below. We also investi- 
gated the dependence on the viscosity value chosen 
for the hot gas. We found that the evaporation ef- 
ficiency is higher for larger viscosity, that means 
the mass flow in the corona is higher and conse- 
quently the mass accretion rate in the underlying 
disk is lower. Our new results are of interest with 
regard to the observed spectra and the spectral 
variations of AGN and galactic black hole candi- 
dates. In Sect. 2 we briefly describe the model. 
The detailed computational results are presented 
in Sect. 3 . In Sect. 4 we discuss the dominant ac- 
cretion form at different accretion rates and the 
possible coexistence of disk and corona in the in- 
nermost region of AGN. Our conclusions are given 
in Sect. 5. 
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2. The model 

We consider a hot corona above a geometrically 
thin standard disk around a central black hole. 
Heat released by friction in the corona flows down 
into lower cooler and denser corona. There it is 
radiated away if the density is sufficiently high. If 
the density is too low, cool matter is heated up 
and joins the coronal gas. We call this process 
evaporation. This goes on until an equilibrium 
density is established. The gas evaporating into 
the corona still retains angular momentum, it will 
differentially rotate around the central object sim- 
ilar to a standard accretion disk. By friction the 
gas looses angular momentum and drifts inward 
thus continuously drains mass from the corona to- 
wards the central object. This is compensated by 
a steady mass evaporation flow from the underly- 
ing disk. The process is driven by the gravitational 
potential energy released by friction in the form of 
heat in the corona. Therefore, mass is accreted to 
the central object partially through corona (evap- 
orated part) and partially through disk (the left 
part of the supplying mass). 

2.1. The equations 

In our earlier study (Meyer et al. 2000a) the 
equations for the corona above the disk are de- 
rived. The approximations used to describe the 
complex physics are discussed there in detail. For 
the application to coronae near the black hole we 
now take additionally into account different ion 
and electron temperature and hence separate en- 
ergy equations for ions and electrons. Also Comp- 
ton cooling is considered. 

We list the equations as follows. 

Equation of state 



P = TliKTi + HeKTe 



|m + Te), (1) 



where we assumed a standard chemical composi- 
tion {X = 0.75, y = 0.25) for the corona. The 
average molecular weight is then /x = 0.62. For 
simplicity we take the number density of ions rij 
equal to that of electrons rie, which is strictly true 
only for a pure hydrogen plasma. 
Equation of continuity 



dz 



{pVz) = -pVr- 



2z 



H -|- z- 



:PVz 



(2) 



^-component of momentum equation 

dvz _ dP GMz 
dz 



dz 



(r2 + ^2)3/2 ■ 



(3) 



For the derivation of the above two equations 
see Eqs. (39) and (34) in Meyer et al. (2000a). 
The term — ^^^^2 P^z is due to the gradual change 
of the channel (1 + z^/r^ the channel cross sec- 
tion) of the ascending gas from a cylindrical to a 
spherical expanding shape at large height z. 

There are now two energy equations (compare 
Eq. (46) Meyer et al. (2000a) for the one- 
temperature case), one for ions and one for elec- 
trons, and the cooling and heating processes in 
such a hot corona have to be included. 

1) For ions the energy balance is between vis- 
cous heating, cooling by collision with electrons 
and radial and vertical advection. 



dz 



PiVz 



<_Pi. _ GM A 
-I Pi (r2-|-z2)V2y 



— PiVr 

2z 
r^-\-z^ 



7 Pi 
7-1 Pi 



GM 



(r24-z2)i/2 
„ ("f ^ I 7 A GM 

"^[2 ^ 7-1 Pi (r-2 + 22)l/2 



(4) 



where the radial velocity is approximated as Vr ~ 
—aVg/V,r {Vs sound speed). The friction is as- 
sumed to be proportional to the total pressure, 
and thus the viscous heating rate per unit vol- 
ume is q-vis = ^aPfl as appears in the equa- 
tion. Here we assume that all viscous heating 
goes to ions as usual (though there is a contro- 
versy on possible viscous heating on electrons (e.g. 
Bisnovatyi-Kogan & Lovelace 1997; Quataert & 
Gruzinov 1999)) and the ion energy is transferred 
to electrons by Coulomb coupling alone. Thus, 
the exchange rate of energy between electrons and 
ions for non-relativistic approximation is (Stepney 
1983), 



lie = (f )'^' If" In AaTcnen,(«Ti - kT,) 



with 



3.59 X lO-^^netiiiTi - Te)^^f^, 



n = — ^ i + — 77^;, 



S72~ 



(5) 



and G gravitational constant, nip and rrie proton 
and electron mass, Thomson scattering cross 
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section, c speed of light, k Boltzmann constant, 
and Coulomb logarithm In A taken as 20. 

2) For electrons the energy balance is deter- 
mined by the processes of heating by collision with 
ions, cooling by Bremsstrahlung, Compton cooling 
and vertical thermal conduction (the radial one is 
neglected since dT^jdr <C dTe/dz), as well as ad- 
vection. 



7 Pe 

7-1 Pc 



GM 



(r2+z2) 
Qie ^e^i-^(^e) ^cornp 

-r^feVr 7-1 (r 



T72) 



2 -,1/2 



GM 



2 -,1/2 



(6) 

where neniL{Te) is the bremsstrahlung cooling 
rate and gcomp the Compton cooling rate. 



4KTe 

2 



(7) 



with Teff the effective temperature of the main 
body of the disk, 



iGMMd 



here we have added a factor 1/2 to the the energy 
density of a black-body photon field in Eq.(7) since 
photons from the underlying disks only cover half 
of the sky of electrons in the optically thin corona 
(compare with Rybicki & Lightman 1979), and Fc 
the thermal conduction (Eq. (49) in Meyer et al. 
(2000a)), 



dz ' 



(8) 



with a radiative constant, a Stefan-Boltzmann 
constant, = 10~^g cm s~^K~^/^ for fully ion- 
ized plasma (Spitzer 1962), and Mj. mass flow rate 
in the disk. 

Prom Eq.(4) and Eq.(6) we obtain an en- 
ergy equation in the same form as for our one- 
temperature model, but with a new term concern- 
ing the Compton cooling. 



dz 



7 p 

7-1 p 



GM 
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= faPfi - n,Ji,L{T, 

PVZ + —J - (^2+^2)1/2 ) + 



2z 
r'-^+z'-^ 



F, 



(9) 



For convenience the two energy equation Eq.(4) 
and Eq.(9) are used in the numerical calculations 
instead of Eq.(4) and Eq.(6). 

We want to point out that the one-zone approx- 
imation is taken for deriving the above equations. 
For a given distance r we choose a region between 
ri and r2 so that 



7r(r^ 



and 



n r' ■ 27rr'dr' 



Ti = 0.72r and r2 



(10) 



(11) 



I.e. ri = U.i'Zr and r2 = 1.23r. In the one- 
zone region the radial advection is taken related to 
the difference between mass inflow from the outer 
neighboring regions and outflow towards the inner 
neighboring regions, 



{2TrrAzpVr)r2 ~ {'2TrrAzpVr)r 
27rrAz(r2 — ri) 



pVr 



(r2 — n) r 

(12) 

This approximation for the net radial flow requires 
the mass outflow larger than the inflow. In the 
following sections we will see that the evaporation 
rate reaches a maximum and doesn't increase any 
more (see Fig. 3) inside the maximal evaporation 
region. This implies our results are inconsistent 
in the inner region and we need to reconsider the 
radial flow. For this we study the radial flow more 
accurately in Sect.4.1. It turns out that the effect 
of the radial advective loss is very small to the 
inner structures, the corresponding terms in the 
continuity and energy equations are not very im- 
portant. Therefore, our results based on one-zone 
approximation are still consistent. 

2.2. Boundary conditions 

The 5 differential equations, Eqs.(2), (3), (4), 
(8) and (9), contain 5 variables P{z), Ti{z), Te{z), 
m{z){= pvz) and can be solved with 5 boundary 
conditions. 

At the lower boundary zq (the interface of disk 
and corona) , there is no heat flux and the tempera- 
ture of the gas should be the effective temperature 
of the accretion disk. Former investigations (Liu, 
Meyer, & Meyer-Hofmeister 1995) show that the 
coronal temperature increases from effective tem- 
perature to 10^'^K in a very thin layer and thus 



2 
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the lower boundary conditions can be reasonably 
approximated (Meyer et al. 2000a) as, 

Te = 10^-^ K, Ti = Te, and = -2.73x lO^P at z = 

(13) 

At infinity, there is no artificial confinement and 
hence no pressure. This requires sound transition 
at some height z = zi, which is taken as our upper 
boundary. There is no heat flux from/to infinity 
either. We then constrain the upper boundary as, 

= and = = p/p = (j^+rj &t z = zi 

(14) 

With such boundary conditions, we have to as- 
sume lower boundary values for P and rh to begin 
the integration along z, if the trial P and rh fulfill 
the upper boundary conditions, we find the solu- 
tions. 

3. Computational results 

We solved the difi'erential Eqs.(2), (3), (4), (8) 
and (9) using the Runge-Kutta method. The 
technical details of the integration procedure are 
similar to those in our previous study (Meyer et 
al. 2000a). In the present investigation standard 
parameters were taken for the black hole mass 
M = IO^Mq and the viscosity coefficient a=0.3. 
The Compton cooling by photons depends on the 
radiation from the underlying disk and hence on 
the mass flow rate in the disk. This is an ad- 
ditional parameter on which the evaporation ef- 
ficiency depends and this parameter is different 
in different astrophysical objects. In order to 
study the coronal properties without involving the 
mass accretion rate in the disk, we first investigate 
the influence of viscosity without Comptonization, 
that means we assumed a very low mass accre- 
tion rate in the disk. The results are presented 
in Sects. 3.1 — 3.5. The computational results for 
the effect of Compton scattering are discussed in 
Sect. 3.6, the relevance for corona around stellar- 
mass black holes in Sect. 3.7. 

3.1. Vertical structure of the corona 

Fig. 1 shows the vertical structure of the corona 
at distances r = 3400rg and r = 340rg (lO^^cm for 
M = IO^Mq), where rg = 2.95xlO"(M/lOSM0)cm 
is the Schwarzschild radius, andT^ir = GMrUpfi/Kr = 
3.37 X 10^'^{r/rg)-'^K the virial temperature. The 



variation with height z of the variables, P{z), 
rh{z), Fc{z), Te{z), Ti{z), and Vz/Vg, is similar 
to that in our previous computations. The tem- 
perature, pressure, density keep almost invariable 
in the corona form z ~ r to the upper boundary 
{z < 3r). Steep changes occur in the lower corona 
(z < r), where the temperature increases from a 
chromospheric value (Tg =Ti » 10^-^K) to a coro- 
nal value (Tj ~ (0.2 — 0.3)Tvir) and the pressure 
decreases rapidly. The density, p oc P/T, drops 
even more steeply than temperature in the lower 
corona. 

At large distance r ~ 3400rg electrons in the 
corona are well coupled with ions, and thus elec- 
tron and ion temperatures are nearly the same, 
Te « lO^K. At a distance close to the black hole, 
Coulomb coupling between ions and electrons be- 
comes weak due to low density and high tempera- 
ture, and the temperature in ions (e.g. Ti w lO^^K 
at r = lOrg) is about 100 time higher than that in 
electrons. 

3.2. Distribution of the temperature ver- 
sus distance 

Fig. 2 shows the change of the coronal tem- 
peratures with distance. The ion temperature in- 
creases towards the central black hole, from val- 
ues of about ^T^ir in the outer region (e.g. Tj « 
1.8 X lO^K at r = 3400rg) to iTvir close to the 
black hole (e.g. Ti « 10"K at r = lOr^), while the 
electron temperature hardly increases beyond the 
value of Te ^ lO^K from the distance r ^ lOOr^ 
inward to the black hole. The temperature distri- 
bution of ions differs from that of electrons since 
the viscous heat of ions can hardly be transferred 
by collisions to electrons and is only transferred to 
the entropy of ions. The fraction of heat turned 
into entropy becomes larger at small distance. The 
electrons, however, can efficiently cool by ther- 
mal conduction and radiation and hence keep their 
temperatures almost constant at decreasing dis- 
tance in the inner region. Such distributions of 
ion and electron temperature versus distance in 
the corona are almost the same as those in a typ- 
ical ADAF (Narayan & Yi 1995; Nakamura et al. 
1997; Manmoto, Mineshige, & Kusunose 1997). 
The evaporation of mass is inefficient in the inner 
accretion region and the viscous heat is transferred 
to entropy. The corona is essentially identical with 
an ADAF. 
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For comparison the temperature distribution 
resulting from our previous one- temperature (1-T) 
model is also plotted in Fig. 2 (dotted line). Prom 
the outer regions to the distance of most efficient 
evaporation (r ~ 340r£,) the 1-T temperature lies 
between the temperatures of ions and electrons of 
two-temperature (2-T) model, actually the three 
temperatures are not much different. In the inner- 
most region the one-temperature solutions are not 
adequate anymore unless there exists much more 
efficient coupling than Coulomb coupling. 

3.3. Distribution of evaporation rate ver- 
sus distance 

We define the mass evaporation rate from our 
one-zone region as Mevap(^) = 27rr^mo(r), where 
mo(r) is the evaporation rate from imit surface 
area at the interface of the corona and the chro- 
mosphere (i.e. the lower boundary). Fig. 3 
shows the rate of mass evaporating from the disk 
to the corona versus distance. The evaporation 
rate Mevap('') increases with decreasing distance r, 
reaches a maximum (hump) of A^hp = 0.012AiEdd 
(MEdd = lOiEdd/c^ Eddington rate) at Vhp ~ 
340rg, and decreases towards the black hole. The 
occurrence of evaporation hump is an important 
characteristic of our disk corona model. The max- 
imum of the coronal flow rate is caused by the 
change in the physical processes that remove the 
frictional heat. At large distances the coronal heat 
is balanced by the inward advection and wind loss. 
Downward heat conduction and subsequent radi- 
ation in the denser lower region play a minor role. 
Thus, the evaporation is weak. At small distances 
thermal conduction removes an increasing part of 
the viscous heat and becomes dominant with rising 
temperature. This downward conducted viscous 
heat then evaporates the cool gas at the chromo- 
sphere and leads to efficient radiation. Therefore, 
the mass evaporation rate in the corona increases 
continuously with decreasing radius until the ra- 
dius Thp is reached. At distances r < rhp, the 
electrons are no longer well coupled with ions in 
the corona and downward thermal conduction is 
less important; on the other hand, since the den- 
sity increases towards the black hole, the heat at 
the lower coronal layer can easily be released by 
radiation in the innermost region, and thus evap- 
oration is not active. 

Also shown in Fig. 3 (dotted line) is the dis- 



tribution of the evaporation rate Movap('') from 
our previous 1-T model (a = 0.3). The maximum 
lies at the same distance from the black hole, but 
its value is one time larger than that of the 2-T 
model. The difference of evaporation efficiency be- 
tween the 2-T and the 1-T model originates from 
the temperature difference, the slightly lower elec- 
tron temperature in the 2-T model leads to less 
thermal conduction and resultant evaporation. A 
small part of heat remains in ions to keep the ion 
temperature slightly higher than the temperature 
in the 1-T model instead of going to electrons and 
evaporating disk gas. 

The gas evaporated from the disk into the 
corona flows inward through the corona towards 
the black hole, except for a small fraction leav- 
ing by wind. Therefore the coronal mass flow rate 
at each distance r is the amount of gas evapo- 
rated from the disk at distances further outward 
up to r. Since the mass evaporation rate has a 
steep hump, the coronal flow at radius r outside 
the hump is approximately A'/ovap(''), while at ra- 
dius inside the hump it is a little bit higher than 
the maximal evaporation rate Mhp- 

3.4. Energy flows in the corona 

For the local energy balance we consider the 
processes of viscous heating, radiative cooling, 
radial advection and vertical advection (energy 
brought into the corona by gas evaporated from 
the disk and removed in the wind). Considering 
the coronal structure at a given distance r from the 
black hole we integrate the energy loss by radia- 
tive cooling and the loss by vertical and sidewise 
advective outflow over z, the height of the corona. 
It is an interesting feature how these losses, nor- 
malized by the release of frictional energy, vary 
with distance r. We plotted the values in Fig. 4. 
It shows that the fraction of viscous heat which is 
lost by radiation is larger at small distances than 
at large distances, while advection is dominant in 
outer the region (r ~ 3400r(,). It takes away the 
major fraction of viscous heat there. At the dis- 
tance r ^ 340rg, where evaporation rate reaches 
its maximum, advection and radiation both re- 
move half of the viscous heat respectively. In the 
inner region the density in the layers at the bottom 
of the corona is large, thus the radiation becomes 
efficient and takes away a larger fraction of viscous 
heat than advection. This feature seems to be in 
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contrast to the typical ADAFs, in which the flows 
are radiation dominated further out and become 
more and more advection dominated with decreas- 
ing distance until finally the advcctivc cooling bal- 
ances the viscous heating in the innermost region 
(Nakamura et al. 1997; Manmoto et al. 1997). 
But with a corona above the thin disk evaporation 
leads to mass exchange and therefore the flow is 
not a constant radial mass flow in the corona as 
in an ADAF. 

The detailed calculations of the vertical struc- 
ture allow to study the difference in the energy 
flow in the upper corona as compared to that in 
the lower layers (compare also Fig. 4 in Meyer 
et al. 2000a). Radiation and thermal conduc- 
tion dominate the energy processes in the very low 
layer, the heat conducted downwards by electrons 
and subsequently radiated away in the lower lay- 
ers is provided by the viscous heat produced in the 
whole region from the chromosphere to a height of 
2; ~ r. In the upper layers the pressure is very low 
and the small amount of viscous heat produced 
there is almost completely advected. Thus, the to- 
tal radiation (vertically integrated over the whole 
corona) predominantly originates from the lower 
coronal layers. Considering the layers from height 
2: = r to the upper boundary, the viscous heat is 
advected, as shown in Fig. 4 (dotted lines). The 
very small fraction of radiative loss (e.g. a few 
thousandths at r ~ 340Grg), decreases towards the 
black hole. In this sense, the corona is advection- 
dominated, consistent with the ADAF except that 
our corona model also includes vertical advection 
flow. 

Comparing the energy flow in OTir computed 
corona with an ADAF, we find that the difference 
between the disk corona and the ADAF is in the 
supply of accretion gas. In the ADAF, accreting 
matter streams inward from the further outward 
adjacent ADAF region, where the temperature is 
also high. While for the coronal flow above the 
disk, at r > r^p, the accreting mass is mainly sup- 
plied from the underlying disk with a temperature 
far lower than that in the corona and the evapo- 
rating gas is heated to a high temperature before 
it accretes inwards. This needs to increase the en- 
tropy and the potential of the gas by "absorbing" 
viscous heat, and meanwhile the heat is partially 
radiated in the dense lower corona. If we consider 
the whole vertical region, the lower and upper 



coronal layers, the average temperature should be 
lower than that of an ADAF with the same mass 
supply rate, and the average radiative cooling is 
more efficient due to the structure of the lower 
coronal layers. After the essential increase of tem- 
perature and corresponding changes of the other 
quantities in the lower corona the quantities are 
nearly constant in the upper layers, where the gas 
density is very low and the temperature is rather 
high. Such a corona is radiation- inefficient, almost 
all of the viscous heat is advected with the accre- 
tion flows (dotted lines in Fig. 4). These typical 
coronal layers are actually like an ADAF, which 
can also be seen from the temperature distribu- 
tion shown in Fig. 2. 

3.5. Effect of viscosity 

The dependence of the evaporation rate on the 
viscosity is shown in Fig. 3. We have chosen a 
relatively large viscous coefficient, a = 0.9, for 
comparing with the results for our standard value 
a = 0.3. As already predicted by an earlier analyt- 
ical study of this dependence, the influence on the 
evaporation efficiency is strong. For a larger vis- 
cosity in the corona the evaporation rate is higher. 
The radial location of the maximal rate is shifted 
inwards for higher viscosity values. The maximal 
rate for a = 0.9 is 0.14MEdd) about 10 times larger 
than that for a = 0.3, and the maximum occurs 
at a distance about 30 Schwarzschild radius, about 
10 times smaller than that for a = 0.3. The rea- 
son is that the larger viscosity leads to more ef- 
ficient heating and an enhanced mass flow in the 
corona. The increased viscous heat is then con- 
ducted downward by electrons and has to be ra- 
diated away by more efficient evaporation of the 
gas in the disk atmosphere. On the other hand, 
faster accretion in the corona tends to decrease 
faster the number density of the corona and re- 
quires more mass supply by evaporation to keep 
the energy balance at the interface of corona and 
disk. Therefore, the increase in viscosity largely 
increases mass evaporation. Similar results were 
obtained in an investigation using the 1-T model 
(Meyer-Hofmeister & Meyer 2001), where small 
viscosity is assumed. 

As an important consequence of this depen- 
dence on viscosity, the distance where the mode 
of accretion via the thin disk changes to an ADAF 
is also affected. For a given mass inflow from the 
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outer part the disk can extend farther inward for 
smaller values of viscosity. Therefore, the emerg- 
ing spectra should be different if two systems have 
a similar mass accretion rate but different viscos- 
ity in the hot coronal gas (for details see Meyer- 
Hofmeister & Meyer 2001). But it is difficult to 
draw detailed conclusions from the observations. 

3.6. Effect of Compton cooling 

At small distance from the black hole Compton 

cooling affects the structure of the corona. The 
Compton cooling term is included in the energy 
equation (9). Considering the processes which 
enter in the energy balance, Compton cooling is 
the only process which causes a dependence of 
the coronal structure on the accretion luminos- 
ity and therefore the mass flow in the underlying 
disk. Otherwise the equilibrium between corona 
and disk relies only on the existence of cooler gas 
below. The dependence on the mass flow in the 
disk adds one more parameter to the evaluation 
of coronal structure and the evaporation rate. To 
eliminate this problem we derived the results as- 
suming a negligible Compton cooling (correspond- 
ing to very low mass flow rate in the disk) in the 
foregoing sections . 

To investigate the effect of Compton cooling 
in a disk around a supermassive black hole with 
a mass of IO^Mq we assumed a rate of Mjn = 
O.O5M0yr~^ (« 0.02MEdd), corresponding to a 
bolometric luminosity L = 2.8 x lO^^ergs s~^ (as- 
suming an energy-conversion coefficient of 0.1). A 
part of the mass flow goes to the corona through 
evaporation. So at a given distance the net mass 
flow rate in the disk should be 

Md{r) = Min-M,{r), (15) 

where Mc(r) is the integrated evaporation rate 
from outer edge to the distance r. This relation al- 
lows to evaluate the energy density of soft photons 
from the underlying disk to be scattered in the 
corona by assuming blackbody radiation. We cal- 
culated the coronal structure and the evaporation 
rate for different distances from large distance to- 
wards the central black hole. In Fig. 5 these results 
for the assumed total mass flow rate 0.02MEdd 
and O.OlMEdd are shown together with the results 
without taking the Compton cooling into account 
(corresponding to a very low mass flow rate). As 



expected the Comptonization has a strong influ- 
ence on the structure of the corona if the mass flow 
rate in the disk is high. Since the mass flow rate 
at distance r in the disk depends on the evapora- 
tion in the further outward regions, whether the 
Compton cooling affects the evaporation rate at 
r depends on the accretion in the outward disk 
plus corona. At large distances, e.g. r ~ lO^r^, 
the corona is hardly affected by Compton scatter- 
ing, while in the inner region, Compton cooling is 
so strong that it overwhelms Bremsstrahlung and 
removes viscous heat efficiently, there is no much 
need for vertically advcctive cooling and hence lit- 
tle gas evaporates from disk to corona. For the ex- 
ample with mass accretion rate 0.02MEdd shown 
in this figure, the evaporation rate at r = 34Grg 
is half of the value without Comptonization, and 
drops very steeply inwards, reaching a value of 
^ 15 times less than the non-Compton one. In 
other words, strong Comptonization in the inner 
region restrains the evaporation. The corona then 
becomes rather weak. 

An interesting feature caused by Compton cool- 
ing is the possible shift of the location of the max- 
imal evaporation rate to a larger distance from 
the black hole. Computations for a series of mass 
flow rates show that from a larger total mass flow 
rate 0.02MEdd to those for negligible mass flow 
in the disk, the evaporation rate at given r con- 
tinuously increases and location of the maximal 
evaporation moves inwards, as shown in Fig. 5. 
For a total mass flow rate from the outer regions 
of Min = 5.0 X 10~^MEdd, the Compton cool- 
ing does not affect the evaporation. For a rate 
Min = O.OlMEdd) there is only very little influ- 
ence by Compton cooling (see the dotted line in 
Fig. 5), at the distance r = 340rc, where the evapo- 
ration rate is maximal, there is already no enough 
mass flow to supply the efficient evaporation. For 
a rate 0.02MEdd the changes are correspondingly 
larger, and the evaporation rate at r = 340rg is 
small and already inside the location of the max- 
imal evaporation rate which moved outward to a 
larger distance. In our previous papers (Liu et 
al. 1999, Meyer-Hofmeister & Meyer 1999) we 
suggested that the transition from the thin disk 
accretion to a coronal flow/ AD AF is caused by 
the mass evaporation. The new results includ- 
ing Compton cooling by the disk radiation which 
show a decrease of the evaporation efficiency in 
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the inner disk regions imply the transition in the 
accretion mode for a given mass flow rate in the 
disk (a rate below the maximal evaporation effi- 
ciency) at a smaller distance from the black hole. 
For higher rates the change between the hard and 
soft spectral states is expected, which occurs at a 
lower mass inflow rate than that without Compton 
cooling. 

The effect of Comptonization obviously de- 
pends on the total mass accretion rate. Since large 
accretion ffow in disk provides strong seed photons 
to be Compton scattered in the corona, the con- 
ductive heat at the lower layers of the corona can 
be appropriately released by Compton radiation, 
the evaporation will decrease or stop from some 
distance inwards. At very large accretion rates, 
our computations show that the Compton cool- 
ing is so efficient that either sound transition at 
the upper boundary is inaccessible, or the corona 
needs artiflcial heat from the upper boundary, in 
other words, no evaporation solutions are found. 
This indicates that probably part of the coronal 
gas falls back to the disk by efScient Compton 
cooling in the innermost region. 

3.7. Corona around stellar-mass black 
hole 

Since the physics of a disk corona around a 
stellar-mass black hole is same as that around a 
galactic black hole, we performed with the same 
code numerical calculations for accretion onto a 
stellar-mass black hole of M = QMq to study the 
influence of evaporation for adequate mass flow 
rates. We obtained similar results for both vertical 
and radial structure. If the distance is measured 
in Schwarzschild radius rg and the mass accretion 
rate in Eddington rate, the radial distribution of 
the evaporation rate and the location of the maxi- 
mum rate are the same as those for a supermassive 
black hole. The values of temperature at any dis- 
tance (in Schwarzschild radius) for both systems 
are also the same. As for the vertical structure, the 
quantities are mass independent as long as we scale 
the size (distance r and height z) in Schwarzschild 
radius and multiply the quantities F^, P, p, m 
by the mass of the black hole. This is essential 
since our equations and boundary conditions are 
invariant against scale transformations that leave 
temperature and velocities constant. Therefore, 
the coronal structures are scale free with respect 



to the mass of central black holes. 

Dimensional analysis of the Compton cooling 
term and its related equations show that the coro- 
nal structures are still independent on the mass 
of the black hole even if the Compton cooling is 
included. This is easily understood by compar- 
ing the Compton cooling term and viscous heating 
term in the energy equation. If we scale the dis- 
tance in Schwarzschild radius and the disk accre- 
tion rate in Eddington rate we obtain the depen- 
dences of Qyis and (?comp on the black hole mass as, 
gvis oc Pfi oc PM-\ gcomp oc j^PM-\ Ob- 
viously the Compton cooling term has the same 
dependence on M as the viscous term, implying 
that the results are also mass independent when 
Compton cooling is included. If we take the rate 
O.OlMEdd (O.O25M0yr-i for a IO^Mq black hole) 
as upper limit for which the corona is not af- 
fected by Compton cooling (the true value should 
be slightly smaller as shown in Fig. 5), we de- 
rive an accretion rate of 2.5 x lO~^'^M0yr~^ (or 
9.5 X lO^'^g s"^) for a 6Mq black hole for the 
limit. This rate is a little less than that expected 
for X-ray transients at spectral transition. This 
means the hard/soft spectral transition of X-ray 
transients is hardly affected by Compton cooling. 
The luminosity for which spectral transitions oc- 
cur would be only slightly lower with Compton 
cooling included. 

Despite the similarity between the corona above 
a stellar mass and that above a supermassive black 
hole (quantities scaled as described) the spectra 
would be different even if the mass accretion rates 
in unit of Eddington rate are the same. The rea- 
son are the different temperatures in the cool disk 
around the black holes. The situation is similar to 
that of the inner ADAFs with outer thin disks. 

4. Discussion 

4.1. Validity of the approximation of ra- 
dial advective loss 

In Sect. 2.1 we already pointed out the radial 

advection problem. For our computation of the 
coronal structure in a given zone of radial extent 
r2 — ri ss 0.5r, we assumed that the mass en- 
tering from outer adjacent region is less than the 
mass evaporating into this zone, and thus we ap- 
proximated the radial mass loss ^^^r^ ~ —^pvr- 
Similar approximation was taken for the energy 



9 



advcction. Computations based on this one-zone 
model show that evaporation does not increase any 
more in the inner region. This means that our ap- 
proximation is justified only for distances from the 
black hole larger than rhp- At smaller distances 
the mass flow in the corona should also include the 
mass entering from outer adjacent region besides 
the evaporation flow, and the radial advective loss 
has to be taken consistently with the mass flow 
change in the corona. 

At a coronal region inside the hump, an inflow 
of ^ Mevap(?'hp) supcrimposes on the evaporation 
flow, increases the coronal density and pressure. 
As a result, the gradient of radial mass Hux—^pVr, 
which is contributed by the evaporation mass in 
the one-zone region, should be formally decreased 
by a factor of ri{r) to compensate the increase of 
the density, i.e.. 



1 drpVr 
r dr 



--r}{r)pvr 



with 



?7(r) 



Mevap(r) 
Mc(r2) -f McvapW 



(16) 



(17) 



where Mc{r2) is the mass flow coming from the 
outer edge of one-zone region and Mevap('') the 
mass flow evaporating into the one-zone region. 
The net energy flow is also multiplied by r?(r). 

Computations have shown that the coronal flow 
comes mainly from the dominant evaporation re- 
gion, i.e. the hump zone. As a rough estimate, we 
simplify ?7(r) to 



/ \ -^^cvap 

7?(r) ~ 



(r) 



Mc(rhp) 



(18) 



Replacing Mevap(?') in Eq.(18) by 2'Kr'^rno{r), we 
recalculate the coronal structure with the new ra- 
dial mass flux and energy flirx, the evaporation 
rate along distance is plotted in Fig. 6. Prom the 
flgure we see that the evaporation by including 
the inflow from outer neighboring region is still 
inefficient in the inner region. The characteris- 
tic feature of hump-evaporation efficiency along 
distances doesn't change. Since is very small 
at small distances far from the hump region, the 
dependence of the evaporation rate on the radial 
advective loss is quite weak. The coronal flow 
inside the hump is the same order of magnitude 
as the hump flow. The evaporation rate and the 



coronal flow seem to depend more strongly on 
the detailed cooling and heating processes, such 
as Comptonization, additional heating (large a). 
On this more consistent global investigations arc 
needed but they are beyond the scope of this pa- 
per. 

4.2. Accretion via corona, disk, or both? 

We investigated the hot corona lying above a 
cool accretion disk and the mass and energy ex- 
change between the corona and the disk. Con- 
tinuously mass evaporates from the disk into the 
corona and flows in the corona towards the cen- 
tral black hole. At large distance from the black 
hole only accretion via a thin disk occurs. Further 
inward (at distances of 10^ to 10^ Schwarzschild 
radii for X-ray transients in quiescence or low lu- 
minosity galactic nuclei) the accretion can change 
to the mode of an ADAF (e.g. Narayan, McClin- 
tock, & Yi 1996). This transition occurs where 
the evaporation rate exceeds the mass inflow rate 
in the cool disk. Prom that distance inward the 
thin disk is depleted and the coronal flow fills the 
inward region up to the last stable orbit. Conse- 
quently, the accretion changes from thin disk out- 
side to coronal flows/ ADAFs inside. The transi- 
tion radius is determined by the values of mass 
evaporation rate and mass flow rate in the disk. 
Actually outside the transition radius both exist 
together, the cool disk and the corona above. Such 
a truncation of the thin disk only occurs if the 
mass supply rate in the disk is less than the max- 
imal evaporation rate. If the mass inflow rate is 
higher, the disk cannot be fully depleted and we 
expect an interior disk below the corona inside the 
distance of maximal evaporation. Summarizing 
these interactions, we claim that the mass supply 
rate Mjn, compared to Mcvap, determines whether 
accretion via the corona or the disk, or both dom- 
inate: 1) At low accretion rate Mi„, a large part 
of the disk is depleted by mass evaporation and 
flUed by coronal accretion flows/ ADAFs, so accre- 
tion is dominated by hot coronal flows/ ADAFs. 
2) At a rate around the maximal evaporation rate 
the disk truncation reaches its innermost possible 
distance, i.e. the hump distance, and an inner 
corona/ ADAF coexists with an outer disk. 3) At 
a higher rate, the disk extends to the last stable 
orbit, and the corona, which is negligible in the far 
outward parts, could be also weak at innermost re- 
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gion due to the efRcient Compton cooling. In such 
a way, the thin disk exists at all distances with a 
corona lying above. 

Up to now the evaporation model was applied 
to the situation of a thin disk truncated at a cer- 
tain distance, that is for mass flow rates Min less 
than the maximal A'/cvap- (Meyer ct al. 2000b, Liu 
& Meyer-Hofmeister 2001). Now the new compu- 
tation of coronal structure for inner regions and 
the evaporation efficiency there connect to the in- 
teresting questions how the radiation from these 
coronae with disks underneath affects the spec- 
trum. 

We have shown that the evaporation rate 
strongly depends on the viscosity value, especially 
that the maximum rate is changed and this maxi- 
mum rate occurs for the higher viscosity at smaller 
distance than for lower a. The larger the viscosity 
value is, the smaller is the inner region with an 
ADAF. Furthermore, the Compton cooling can 
affect the evaporation. If the mass flow rate Min 
is high, the Compton cooling reduces the evapora- 
tion, and due to the reduced maximal evaporation 
rate the disk extends further inwards and may 
even extend to the last stable orbit. Therefore, 
the dominant accretion mode and the transition 
radius from outer disk accretion to an inner ADAF 
sensitively depend on both, viscosity and mass 
supply rate to the disk. The dependence of the 
transition radius on the mass accretion rate is dif- 
ferent from that in our previous investigation if 
Compton cooling becomes important. 

4.3. Corona + thin disk in AGN? 

In AGN hard X-ray tails and iron lines are 
observed (e.g. Pounds et al. 1990; Matsuoka et 
al. 1990; Tanaka et al. 1995), implying that cold 
and hot gas exist in the vicinity of the central 
black hole. Commonly it is believed that these 
components are corona and disk respectively. Soft 
photons from the disk are Compton up-scattered 
by energetic electrons in the optically thin corona 
and form the hard X-ray power-law spectral com- 
ponent. The iron lines are from the illumination of 
disk matter by X-ray radiations from the corona 
(e.g. Guilbert & Rees 1988; Lightman & White 
1988; George & Fabian 1991). 

The disk corona evaporation model shows that 
such kind of structure can exist self-consistently, if 



the mass accretion rate and the viscosity have ap- 
propriate values. Considering the high luminosity 
in AGN, Comptonization is quite strong and the 
evaporation is restrained. Near the black hole part 
of the coronal gas is expected to cool down and 
fall back to the disk. To keep the corona above 
the thin disk at this innermost region, very largo 
viscosity is required as is indicated in Fig. 3. This 
implies that a large fraction of accretion energy 
needs to be released in the corona. Here wc raise 
the same question as Merloni & Fabian (2001b) 
did, that how the corona stores and releases so 
much energy. 

In luminous AGN (i.e. large M/MEdd)) corona 
can not exist in the innermost region anyhow 
as Compton cooling is so strong that coronal 
plasma completely cools down and falls back to 
the disk. No strong hard X-ray spectral com- 
ponent is expected from corona. This is consis- 
tent with the observational feature that high-m 
objects (e.g., Narrow-Line Seyfert 1 galaxies with 
relatively small black hole mass) show steep hard 
X-ray spectra (e.g. Boiler, Brandt, & Fink 1996; 
Brandt, Mathur, & Elvis 1997), while the rela- 
tively low luminosity AGN (e.g., classical Broad- 
Line Seyfert Is) show hard X-ray tails. 

5. Conclusions 

We studied the properties of disk and corona 

around black holes including the exchange of mass 
and energy between the disk and the corona. The 
decoupling of electron and ion temperatures and 
Compton cooling were taken into account. We 
found that the viscosity parameter a for the coro- 
nal gas has a strong influence on the evapora- 
tion efficiency, for larger a the mass flow in the 
corona increases and the maximum of the evap- 
oration rate is reached at a smaller radius. We 
also found that mass evaporation to the corona is 
largely restrained by Compton cooling at regions 
close to the black hole. Both effects, larger vis- 
cosity and smaller incoming mass accretion rate 
to the disk, lead to the coronal flow/ADAFs ex- 
tending outwards to a larger distance. Coronal 
radiation is weak for high accretion rates, consis- 
tent with the steep hard X-ray spectra observed 
in NLSls and galactic black hole candidates. To 
maintain a corona lying above a thin disk down 
to the last stable orbit in AGN, a large energy 



11 



store/release is required. The results based on 
galactic nuclei are mass independent and therefore 
also applicable to black hole X-ray binaries. 
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Fig. 1. — Vertical structure of the corona at different distances from a central black hole of mass IO^Mq, at 
r ~ 3400rg(left panel) and r = 340rg (right panel). Distributions of the following quantities: Temperatures 
of ions Tj and electrons Tg in units of 0.2Tvir, (TVir virial temperature); pressure P, thermal conductive flux 
Fc and vertical mass flow m = pv^ scaled in their values at the lower boundary. Steep variations of P, Ti, 
Te, and rh occur in the lower layers with z <r, above these transition layers the corona is quasi-constant. 
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Fig. 2. — Temperature distribution versus dis- 
tance from the black hole. Solid lines lower to 
upper: distributions of electron and ion temper- 
ature Te and Tj, and virial temperature Tvir, re- 
spectively. Dotted line: temperature distribution 
of a 1-T model (Liu & Meyer-Hofmeister 2001). 
The ion temperature increases towards the central 
black hole, the electron temperature stays around 
lO^K in the inner region. The distribution result- 
ing from the 2-T disk corona model is almost the 
same as that from a typical ADAF. 
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Fig. 3. — Evaporation rate versus distance. Solid 
lines: results from the 2-T model for a viscous co- 
efficient a = 0.9 (upper line) and a = 0.3 (lower 
line); dotted line: result from the earlier compu- 
tations using a 1-T model, a = 0.3. The mass 
evaporation rate increases more than 10 times 
when the viscosity increases 3 times. The evap- 
oration rate resulting from 2-T model is a lit- 
tle smaller than that from the 1-T model since 
a small part of heat remains in ions to maintain 
a slightly higher temperature and less heat goes 
to electrons (lower electron temperature) and is 
conducted downwards for evaporation. 
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Fig. 4. — Variation of the energy losses with dis- 
tance r from the black hole. Fraction of viscous 
heating going into radiation and fraction removed 
by advcction. Solid lines: rates integrated over the 
total corona from the lower boundary (zo, chromo- 
sphere) to the upper boundary (sound transition). 
Dashed lines: rates integrated over the transition 
layer only, from zq to z = r. Dotted lines: rates 
from z = r to the upper boundary. The com- 
parison of solid and dotted lines show that the 
processes in the lower coronal layers arc respon- 
sible for the higher radiative loss fraction in in- 
ner region and the lower advective loss fraction in 
the outer region. Considering the upper coronal 
layers only, the viscous heat almost flows inwards 
as entropy with little radiation loss (a few thou- 
sandths in outer region), which is quite similar to 
the ADAFs. 
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Fig. 5. — Effect of Compton cooling on the evapo- 
ration. Upper solid line: Evaporation rate without 
Comptonization (corresponds to a very low mass 
flow rate in the thin disk); lower solid line and 
dotted line: including Compton up-scattering of 
soft photons from the disk with a total accretion 
rate of Mn = 0.02MEdd and Mi„ = O.OlMEdd, 
respectively. The gas is partially evaporated into 
a coronal flow and partially accretes through the 
disk. The figure shows that the evaporation can be 
strongly restrained by Compton cooling, depend- 
ing on the mass flow rate in the disk, in particular 
in the inner region. 
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Fig. 6. — Influence of mass inflow from the hump 
region on the evaporation. Solid hne: basic case 
(a = 0.3 without Compton coohng). Dotted hne: 
influence of hump flow. Dashed hne: influence of 
Comptonization. Dash-dotted hne: influence of a. 
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